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ABSTRACT 

The  problem  of  a  thick  cylindrical  antenna  driven  by  a  "delta-function 
generator"  is  investigated.  Numerical  solutions  are  obtained  for  two  dif¬ 
ferent  mathematical  models  —  one  is  driven  by  a  "delta-function  generator" 
on  both  outside  and  inside  surfaces  of  the  antenna,  and  another  is  driven 
only  on  the  outside  surface.  In  both  cases,  current  singularity  near  the 
driven  point  has  been  taken  care  of  and  subsequently  subtracted  out,  before 
the  numerical  solution  is  applied.  For  mtennas  with  circumferences  less 
than  a  free-space  wavelength,  the  results  arc  compared  with  the  experi¬ 
mental  data  obtained  by  Holly  [12].  For  antennas  with  larger  radii,  no 
experimental  data  are  yet  available. 
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I.  INTRODUCTION 


In  the  theory  of  linear  antenna,  the  total  current  distribution  along 
an  antenna  with  axial  symmetry  can  usually  be  obtained  by  solving  an 
integral  equation.  It  is  well  known  that  the  integral  equation  of  the  dipole 
antenna  can  be  written  as  [i] 

J  k  IJz')  K(z-z';  a)  dz'  =  (c  cos  k  z+-^-  VsinkQ|  z|  V  |z|<h 
~h  4o  ° 

(1) 

where 


+7J.  +ik  -%/ (z~z')2  +  (r^  +  a^  -  2racos  0) 

K(z-z';  a)  =  ■—•  f  d0  . .  . - . -  ~  - - - — 

^  \J (z-z1)^  +  (r^  +  a^  -  2racos  0) 


(2) 


kQ  the  free  space  wave  number  and  the  characteristic  impedance 
of  free  space  (i.  e.  ,  120  7 r  ohm).  The  integral  equation  (1),  which  is  sat¬ 
isfied  by  the  total  current  distribution  I^.{z)  on  the  antenna,  corresponds 
exactly  to  the  mathematical  model  of  a  tubular  dipole  antenna  which  has 
a  half-length  h  and  radius  a  (Fig.  la).  Near  the  origin,  z  =  0,  the 
antenna  is  dirven  by  a  so-called  "5 -generator"  of  voltage  V,  such  that 
the  tangential  electric  field  E^(z)  on  both  outer  and  inner  surfaces  of  the 
antenna  is  given  by 

Ez(z)  =  -V6(  z).  (3) 

For  antennas  with  very  thin  radius,  various  methods  have  been 
developed  for  solving  this  integral  equation.  For  example,  King- 
Middleton  in  1946  had  already  obtained  a  second-order  iterated  solu¬ 
tion  [2].  Tai  in  1950  used  a  variational  technique  to  obtain  a  first-order 
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solution  whose  input  impedance  is  stationary  [3],  King  in  1959,  and 
King-Wu  again  in  1964  derived  a  two-term  and  three-term  theory 
for  the  current  [4].  Wu  in  1956  applied  the  Weiner-Hopf  technique  to 
the  problem  of  the  long  antenna  [5].  Duncan-Hinchey  in  1960  used  the 
Fourier-series-expansion  technique  [6]  and  Mei  in  1964  obtained  a 
numerical  solution  for  the  current  distribution  [7].  The  solutions  obtained 
by  these  authors  were,  in  one  way  or  another,  based  upon  an  approxima¬ 
tion  of  the  kernel  K(z-z';  a)  of  the  integral  equation.  These  approxima¬ 
tions  are  satisfactory  only  when  the  radius  of  the  antenna  is  very  small 
compared  with  the  free-space  wavelength  and,  therefore,  cannot  be 
applied  directly  to  the  antenna  with  moderately  large  radius. 

Moreover,  as  pointed  out  recently  by  King  and  Wu  [8],  the  model 
corresponding  to  Eq.  (1)  does  not  agree  precisely  with  what  is  used  in 
the  experiment.  Experimentally,  the  usual  coaxial-line -driven  antenna 
over  a  ground  plane  is  excited  only  from  the  outside  and  not  from  both 
sides  of  the  antenna  (Fig.  lb).  In  order  to  represent  the  experimental 
model,  an  extra  term  has  to  be  added  into  Eq.  (1),  i.  e. , 


J  k^z1)  K(z-z';a)  dz'  =  ( C  cos  k^z  +  y  sin  |kQz  | ^ 

-h 


i27rk^  a 

+  -—V  f  r'dr'K(z-z';r') 
^o 


z'=0’ 


<  h. 


(la) 

2 

As  is  shown  later,  this  extra  term  in  Eq.  (la)  is  in  the  order  of  (kQa) 
and  therefore  can  be  neglected  for  thin  antenna.  Eq.  (la)  then  reduces 
to  Eq.  (1).  This,  however,  may  not  be  true  for  antennas  with  moderately 
large  radii. 
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On  the  other  hand,  with  the  help  of  the  modern  high-speed  computer, 
integral  equations  of  this  type  (i.  e. ,  Eqs.  (1)  and  (la))  can  actually  be 
solved  numerically  by  a  so-called  "approximate  product-integration  tech¬ 
nique,"  even  when  the  exact  kernel,  which  has  a  singularity  at  z  =  z\ 
is  retained.  The  purpose  of  this  report  is  to  apply  this  method  to  the 
thick  antenna. 


II.  DIPOLE  ANTENNA  WITH  2 -SIDED  EXCITATION 


A.  CURRENT  BEHAVIOR  NEAR  THE  DRIVING  POINT 

Let  Z  =  kQz,  H  =  k  h,  A  =  kQa  so  that  Eq.  (1)  can  be  normalized 
as  follows: 

+H 

/  I  (Z1)  K(Z-Z';  A)  dZ1  =  — -  (c  cos  Z+  V  sin  |z  |  );  |z|^H.  (4) 

-H  1  * o  A  y 

The  boundary  condition  at  both  ends  of  the  antenna  requires  the  total  cur¬ 
rent  I^Z')  to  vanish,  i.  e.  , 

It(±H)  =  0. 

In  order  to  solve  Eq.  (4)  numerically,  the  integral  on  the  left  must 
be  approximated  by  a  sum  of  products  of  some  weighting  functions  and 
the  values  of  I t(Z)  at  certain  sample  points.  Hence,  any  singularity  in 
the  current  that  appears  in  the  solution  of  Eq.  (4)  must  be  evaluated  first 
and  subtracted  out  from  the  integral  in  order  to  assure  a  gopd  approxi¬ 
mation.  From  the  analysis  of  thin  antennas,  it  is  well  understood  that 
the  singular  part  of  It(Z)  is  due  directly  to  the  existence  of  the  non¬ 
physical  6 -function  source,  and  that  this  par  ;  of  the  current  is  independ¬ 
ent  of  the  antenna  length.  That  is,  when  H  —  oo,  Eq.  (1)  will  degenerate 
into 

f  I.(Z')  K(Z-Z>;  A)  dZ'  =-^e+llZl  =  M(Z);  -co  <  Z  <  co  (4a) 

^  1  Y 

-00 

Now  define  the  Fourier  transforms  of  It(Z),  M(Z)  and  K(Z-Z';  A)  as 
It(\),  M(\)  and  K(\,A),  respectively. 


2-1 


2-2 


V®  ■  t 

It(\)  e"lXZ  d\ 

(5a) 

-  2i 

V- 

/  °°  It(Z)  e+iXZ  dZ 

-oo 

(5b) 

M(Z)  =  /  MOO  e"iXZ  dX, 

(6a) 

C 

§ 

II 

do) 

C 

f+°°  2*V  e+ilZl  e+iXZ  tfZ 
-00  ^0 

i2V 

1 

(6b) 

^0 

1 

CO 

K(Z,R)  =  J  K(\Jl)  e  "iKZ  d\ 

C 

1  +  00  i  •%  ry 

K(\,R)  /  K(Z,  R)  e  lXZ  dZ 

-00 


(7b) 


Since  the  kernel  K(Z-Z';  R)  is  the  free  space  Green's  function  of  a  ring 
source  which  is  distributed  on  the  surface  of  the  antenna,  the  Fourier 
transform  of  the  Green  function  can  be  substituted  in  Eq.  (7b).  Thus, 


K(\,R)  =  J  ^  H 1  -X2  \/r2  +  A2-2RA  cos  e) 

-7T  / 


de 


,  +00 

=  ^F  /„  dC  £  Hn  ] \R>  J„(B<  V^) 


+in0 


n=-oo 


*THo ’(*>  a/77^) 


(7c) 
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where  (R>  A).  The  branch  cuts  at  X  =  ±1  and  the  path  of  integra- 

< 

tion  C  are  shown  in  Fig.  2.  The  choice  of  these  branch  cuts  is,  in  prin¬ 
ciple,  quite  arbitrary.  One  criterion  is  that  the  solution  obtained  for 
|x|  <  1  on  the  real  axis  should  be  in  agreement  with  the  solution  obtained 
for  |x|  >1,  also  on  the  axis,  after  detouring  around  the  branch  point. 
However,  for  the  convenience  of  numerical  integration,  the  branch  cuts 
actually  chosen  are  the  following: 

(i)  k2  -  1  =  V  \2  -  1 

(ii)  /\J X2  -  1  =  -i|/\/  X2  -  1  j 

(iii)  xj X2  -  1  =  +i|/\/  X2  -  1  j 

(iv)  |  arg  \/  X2  -  l|  <  -~ 


X  >  1  on  the  real  axis; 


on  the  hand  side  of  the 

left  right 

branch  cut; 

on  the  right  hand  side  of  the  light 
branch  cut; 

elsewhere.  (8) 


With  the  branch  cuts  defined  above,  (5),  (6)  can  be  substituted  in  (4a),  and 
use  made  of  the  convolution  theorem.  Then 


2*It(X)K(X,  A)  U 


■‘•J.V'W  „  _  _  ' 

71  0  (1-\2)H^(a'\/  1  -x2)  jq(a  a/  1  -X2) 

I  (Z)=-^-  f  dX  e"iXZ - ; - - ;  - 

c  (1-X2)H^(a^)jo(aV^) 


cos  XZ 

dX - - - - 

(l-x2) H^(aV  1  -x2)  J0(aV  1  -x2) 


(10) 
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Here,  since  the  integrand  behaves  as  — — —  near  the  branch  point  (i.  e. , 

c  in  £ 

+i  0. 

\  =  1  +  e  e  ),  the  integration  around  the  branch  point  vanishes  as  e  -  0. 

Although  the  expression  for  It(Z)  for  all  Z  is  rather  difficult  to 

obtain,  the  leading  term  in  I^(Z)  near  Z  =  0  can  be  obtained  directly  from 

Eq.  (10).  That  is,  if  X  is  assumed  to  be  a  large  number  such  that  \q  »,1, 

A'n/  1  +  »  1,  the  integral  can  be  split  into  two  parts,  one  evaluated 

from  0  -*  X  ,  the  other  from  X.  to  oo.  For  the  integral  from  X.  to  oo,  use 
o  o  o 

can  be  made  of  the  asymptotic  expressions  for  the  Bessel  function. 


*  0 


/  °  <*- 


cos  XZ 


(-iflrA)  / 


( 1  -X 2)  H^1 ) (aV  1  -X2 )  J Ja* 
.....  1 


cos  \Z 


cos  XZ 


_4v_  o  ^ - 

(I-X^H^aV  1-\2)jo(aV  i-x2) 


+  i7TACi(X  Z) 
o 


Here,  Ci(XQZ)  is  the  cosine  integral  whose  behavior  near  \QZ  —  0  is 
known  as 


V  <-Dn 

Ci(XoZ)=  Y+lnX0Z+  2,  - ; 

n=l  2 


(-1 )  (XQZ)‘ 
2n(2n) ! 


The  substitution  of  Eq.  (12)  into  Eq.  (11)  gives 


I  (Z)  ~  in  x  Z  as  Z  —  0. 

t  r  o 


This  logarithmic  term,  of  course,  had  been  observed  early  in  1959  by 
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King  and  Wu  [9].  Also,  since  this  logarithmic  term  represents  a  purely 
local  effect  due  to  the  6 -function,  it  should  be  essentially  unchanged  for 
an  antenna  of  finite  length. 

B,  NUMERICAL  SOLUTION  FOR  THE  CURRENT  DISTRIBUTION  - 
THE  "APPROXIMATE  PRODUCT-INTEGRATION  TECHNIQUE" 


1.  Total  Current  Distribution 


To  solve  the  integral  equation  numerically,  each  half  of  the  antenna 

is  first  divided  into  N  equally  spaced  segments  so  that  each  segment  has 

H 

the  length  2A  where  A  =  ^r-.  Furthermore,  if  Z  =  (n-l)A,  n  =  1,  2, , 
2N+1  (Fig.  3),  the  right-hand  side  of  Eq.  (4)  can  be  written  as  the  sum  of 
N  integrals.  The  nth  integral  in  the  equation  is 


+A 


/  It(Z'+Z2n)[K(Z-Z'-Z2n)  +  K(Z+Z'+Z2n)]  dZ'. 


(14) 


In  this  segment,  let  it  be  assumed  that  It(Z)  can  be  approximated  by 
M 

Vz>  =  I  +  <16) 


s  =  l 


where  6^  is  equal  to  1  if  n  =  1  and  0  if  n  ■+  1.  The  first  part  in  Eq.  (15) 
actually  corresponds  to  the  first  M  terms  in  the  Taylor  series  expansion 
about  the  center  of  this  segment.  Thus,  the  choice  of  M  and  the  accuracy 
of  the  approximation  depends  upon  the  true  solution  I^(Z)  and  its  deriva¬ 
tives.  The  second  part  in  Eq.  (15)  is  a  logarithmic  term  which  exists 
only  in  the  first  segment  and  is  normalized  co  vanish  at  the  end  of  the 
segment. 


;V/WT 

•Wi  I 

ff-i 


fc^l 


¥ 


BLANK  PAGE 
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Now  let  u.  with  j  =  1,  2, ...  M  be  the  M  sample  points  in  the  seg- 

J 

ment  and  let  the  integral  in  (14)  be  expressed  as  a  sum  of  the  products 

of  some  "weighting  functions”  a^)( z)  and  the  values  of  I,(Z)  at  the  sample 

nj  t 

points,  i.  e. , 


/  It(Z'+Z2n)[K (Z-Z'-Z2n)  +  K(Z+Z'+Z2n)]  dZ» 


AT4. 

-  I  4*j(Z)  ¥uj>  +  5nlf<t)(Z)' 

3=1 


where 


i\z)  = 


+A  (Z'+A) 


/  In  — 2X —  [K(Z-Z>-Z  )  +  K(Z+Z'+Z  )]  dZ' 


and  the  a  .  (Z)  are  to  be  determined.  Equation  (15)  is  then  substituted  on 
both  sides  of  Eq.  (16)  to  give: 


I  V^8'1  C  (-ff1  [K(Z.-Z'-Z2n)  +  K(Z+Z'+Z2n)l  dZ' 


-  Y  a(t)(Z)(u.-Z9  )S'H  =  0 

Li  niv  j  2n 

3=1  J 


a  AS“V(t)(Z) 
sn  rns 


-  1  • °:  (18) 
3=1 


fx' l,(Z )  is  called  the  "moment  function."  Since  Eq.  (18)  is  true  for  any 
ns 

asn'  eac^  lerm  summation  is  zero.  Therefore 
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1 

a(t|(Z) 

nl 

(t) 

an2  Z 

YY(Z> 

/  7  \M-1 

C“  i !-) 

• 

• 

• 

_°n> 

• 

* 

• 

Vm'  ' 

(19) 


or  simply. 


*„«;>  ■  p.«> 

where  0^t}  =  {c^  .(Z)}.,  =  {^.(Z)}.,  and  d£>n  =  |(^~ 


u.  -  Z 


2rf 


(19a) 

i-f 


A 


Multiplication  of  both  sides  of  (19a)  by  the  inverse  matrix  of  X  yields 
the  expression  for  c^.  in  terms  of  the  moment  function  jx^.(Z). 

The  error  involved  in  Eq.  (19)  is  essentially  due  to  the  remaining 
terms  in  the  Taylor  series  expansion.  It  is  possible  to  estimate  this 
error  and  to  "throw  it  back"  to  improve  the  solution,  but  this  will  not  be 
discussed  here  [10].  Also,  Eq.  (19)  can  be  simplified  further  if  It(Z)  is 
assumed  to  be  approximated  by  a  lower -order  polynomial  and  M  to  be  a 
set  of  equally  spaced  sample  points.  For  a  parabolic  approximation,  let 
M  =  3  and  m  =  Z2n-1  +  =  Z2n+j-2*  Then 


JA 


Of  _ 

- 


1  1  1 

-1  0  1 

1  0  1 


(20  a) 
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The  substitution  of  (20a)  into  (19a)  leads  to  an  explicit  expression  for 


a  (Z). 

n  y 


^1<Z>  =i[-^2<Z>  +  ■‘"s®] 


_(t)  /lTr\  _  (t)  (t)  (y\ 

n,  2  Z  "  pn,  1(Z)  ■  ^n,  3(Z) 


(t)  fy\ 

Gn,  3  Z 


With  Eqs.  (16)  and  (20b),  Eq.  (4a)  becomes 


N  _3 

1  I^)(z)It(z 

n=l  j=l 


2n-2+j 


(20b) 


•)  +  f(t)(Z)  =  ^  (c  cos  Z  +  —  sin  |z|);  Z  <  H, 


or  simply,  at  Z  =  Z^, 


2N+1 


I  Wm  +  fK  =  (C  ZK+  T  sin  Zk)  • 


m=l 


K  =  1,  2, . . .  2N+1, 


where 


g(t)  _  „(t)  try  \ 

PK1  '  1,  1(ZK' 


S  =  1,  2, . .  .N 


S  =  1,  2, .  .  .N 


p(t)  _  (t)  (7  V 

PK,  2S  "  S,  2  ZK  ’ 

o(t)  _  (t)  /nr  \  ,  (t)  (7  V 

PK,  2S-1  CS,  l'^K' +  S-1,3(ZK'’ 

«(i)  _  Q(t)  \ 

PK,  2N+1  ‘  °N,  3  iv 


Since  the  current  It(Z)  must  vanish  at  the  ends,  ^N+l  can  se^  eclua^  to 
zero.  Equation  (22)  then  involves  only  2N+1  unknowns.  Thus,  if  ^N+l 
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and  p 


-i47T 


K,  2N+1 


cos  Zj£  are  redefined,  Eq.  (22)  can  be  rewritten  as 


follows: 


2N+1 


I 

m=l 


^(t)  r(t) 
rKm  m 


i27T  .  ry  c 

—  V  Sln  ZK  ‘  fK; 


K  =  1,2,...  2N+1  (24a) 


'o 


or  simply. 


(24b) 


where  m,  $t  =  Om>  and  Q  -  V  sin  ZR  -  f^. 

The  multiplication  by  the  inverse  matrix  of  ^  on  both  sides  of  (24b) 
allows  the  values  of  It(Z)  at  the  sample  points  to  be  obtained. 


2.  Outside  and  Inside  Current  Distributions 


For  any  point  (Z,R)  in  air,  the  magnetic  field  due  to  the  current  dis¬ 
tribution  on  the  antenna  can  be  expressed  as  the  curl  of  the  vector  poten¬ 
tial, 

-K  a  +H 

Bq(Z,  R)  =  (VXA)q  =  ~  /  It(Z')  K(Z-Z';  R)  dZ>.  (25) 

o  -H 


Thus,  the  parts  of  the  current,  which  are  distributed,  respectively,  on 
the  outer  and  inner  surface  of  the  antenna,  can  be  determined  from  the 
boundary  condition  on  this  surface, 

I  ,(Z)  =  ±27 tv  aBQ(Z,  R=A±e) 
out  0  0 

in 


_9_ 

9R 


I 


+H 


-H 


It(Z')  K(Z-Z';R)  dZ' 


R=A±e 


as  e  —  0. 


(26) 
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The  difference  between  these  two  currents  is 


IJZ)  =  I  ,(Z)  -  I.  (Z) 
d  out  m 


(27a) 


A  9 

T  81  /H  h{Z[)  K<Z-Z';R)  dZ' 


R=A+e 


+  f  It(Z')K(Z-Z';R)dZ' 

“'ll 


;  as  e  —  0 


R=A-e 


a  +H  r  a  9 

4  f  I.(Z’)  ~K(Z~Z';R)  +WK(Z-Z';R) 

2  -H  1  8R  R=A+6  9R  R=A-€ 


dZ'; 


as  e  -*•  0. 


The  interchange  of  differentiation  and  integration  in  the  last  step  of  Eq.  (27) 

is  not  necessarily  true  for  I  .  or  I.  alone. 

^  out  m 

A  comparison  of  (27)  with  the  left-hand  side  of  Eq.  (4a)  leads  to  the 
definitions 


d" A  /^,vs-l 


=A-e 


4  4  K(Z+Z.+Z2n;R)  R=A+^  +  ^  K(Z+Z.+Z2n;R) 


as  e  —  0 
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f(d)(Z)  = 


i4VA 


+A  (Z'+A) 


f  l.n 


^0 


2A 


A 

2 


^K(Z-Z--Z2n;R) 


R=A+e 


+  _^K(z-z<-z2n;R) 


R=A-e 


A 

2 


^K(Z+Z>+Z2n;R) 


R=A+e 


+  FRK(Z+Z,+Z2n;R) 


R=A-e 


>  dZ»; 


as  6  -  0. 
(29) 


Equation  (27)  can  then  be  transformed  into 


W  * 


e(d)  r(t) 
Km  m 


(30) 


m=l 

where  =  f^(Zv)  and  are  related  to  the  moment  functions  jx^^(Z) 

iv  is.  xvin  ns 

in  the  same  manner  as  in  Eqs.  (20b)  and  (23).  Let  Eq.  (30)  be  combined 
with  Eq.  (27a)  to  obtain 


Xout^ZK^  =  2 
in 


i<t) 


fi(d)  T(t) 
pKmm 


(31) 


C.  EVALUATION  OF  THE  MOMENT  FUNCTIONS 


Define 


H^(S,m)  = 


K(mA-Z';A)  dZ' 


(32) 


2-1? 


^^(s.m)  -- 


/+*  i^rf  1  (-f )  \ik K<^-z') b) 


-  ~  K(mA-Z';R) 


R=A-e 


dZ'; 


R=A+e 
as  €  —  0. 


(33) 


As  a  consequence  of  the  symmetry  of  the  kernel,  i.  e. ,  K(Z,  R)  =  K(-Z,  R), 
the  moment  functions  p^(,s,m)  and  |J.^(s,  m)  have  the  following  proper¬ 
ties: 

(i)  p(s, -m)  =  (-l)s  1  p.(s,m) 


(ii)  |x(l,  0)  =  2p(l,  1) 


(iii)  p(2,  0)  =  0 


A  A 

as 


(34) 


(iv)  |x(3,0)  =4[^(3,  l)-2p(2,  1)+Ht(l,l)]  A 

as  A-~ 


Equation  (34)  holds  true  for  both  subscripts  (t)  and  (d).  A  comparison  of 

(34)  with  the  definition  for  ^(Z)  leads  to 

rij  s 

p(t)'  (d)(Z„)  =  (d)(s,  K-2n)  +  (d)(s,  -K+2-2n).  (35) 

n,  s  xv 

After  the  substitution  of  (35)  in  Eqs.  (20b)  and  (23),  it  is  possible  to  ex¬ 
press  in  terms  of  p^J^(s,m)  explicitly,  i.  e. , 

PK  i  =  Y[^(3,|2-Kl)+H(2-K)^2j2-Kj)+[Ji(3,K)+p(2,K)] 

PK  2S_1  =-|[H(2S-K)tx(2,|2S-K|)+Hi(3j2S-K|)+[x(2,2S+K-2)+M.(3,2S+K-2) 

-  H(2S-K-2)p(2j2S-K-2|)+(x(3,|2S-K-2|)-(x(2,2S+K-4)+fx(3,2S+K-4)] 
PK  2S  =  H(1,|2S-K|)  -  p(3,|2S-Kj)  +  p(l,2S+K-2)  -  ,a(3,2S+K-2) 
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where  H(Z)  is  defined  as  -1  if  Z  <  0  and  1  if  Z  >  0. 

1.  Evaluation  of  |J^(s,m) 

The  substitution  of  Eq.  (7a)  and  Eq.  (7c)  in  (32)  and  the  interchange 
of  the  two  integrations  gives 

P(t)(s,m)  =  ~f  dX  H^aV  1-X2)  Jq(aV  1-X2)  /  AdZ'  *  e+lX(Z,"inA) 

"  C  -A 

=  y/ |Gs(na)  H^(aa/^)  Jq(aV  Hi2)  dX,  (37) 

J  c 

where 

Gs(mA)  =  -i(l-e-i2X^  e-i(m-1,XA;  S  =  1 

=  -^[(ixA-l)+(AA+l)e-i2XA]  e-i(m'1)XA;  S  =  2 

=  ‘(xk)  [^2A2+i2^-2>-(x-2A2-i2XA-2)e'l2XA]e"1*m"1)XA;  S  =  3. 

(38) 

Here,  it  can  easily  be  seen  that  G<^m,  X)  -  -i  as  \  •*  oo,  and  that  Gg(m,-iX) 
is  purely  imaginary. 

In  Eq.  (37),  aside  from  the  logarithmic  singularities  at  the  branch 
points,  the  integrand  has  no  other  singularity.  Furthermore,  as  |xj  —  oo 
in  the  lower  half  plane, 

f  Gg(X)  H^aV  l-x2)  Jo(A\A:r2)  -const e"(m"1)A  l1^  I . 
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Therefore,  the  contour  can  be  deformed  from  C  to  C1  for  m  ^  1, 
H(t)<S,m>  -  i  (f*  +  £  +  jj  +  j£)  ^  Gg(m,X)  j£>  (a  e+  2 


i_ 

2 


Gg(m,\) 


H 


d)i 


(A  e+i7r|A/x^T  ) 


Equation  (39)  is  the  final  expression  for  computation.  The  first  integral 

consists  of  both  real  and  imaginary  parts,  and  is  integrated  within  a  finite 

interval.  The  second  integral  has  only  a  real  part  and  is  integrated  from 

0  —  oo;  the  integrand  decays  at  least  as  as  X.  oo.  In  both  integrals, 

X 

due  to  the  choice  of  the  branch  cut,  the  arguments  of  the  Bessel  functions 


are  all  real. 
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Likewise,  an  expression  for  fj^  can  also  be  obtained  (see  Appendix). 
2.  Evaluation  of  |x^(ni,S) 


The  substitution  of  Eqs.  (7a)  and  (7b)  into  Eq.  (33)  gives 


As  in  Eq.  (37),  the  integrand  in  Eq.  (40)  also  has  two  branch  points  at 
\  =  ±1.  However,  the  integrand  is  finite  at  these  points.  On  the  real 
axis,  as  |\  |  —  oo 


const  — g 


(41) 
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Therefore,  due  to  the  cancellation  of  the  leading  term,  Eq.  (40)  is  inte¬ 
grate  on  the  real  axis  for  all  m  >  1.  However,  this  is  not  true  on  the 
imaginary  axis,  since  as  |x  |  —  oo. 


*  (42) 

The  (+)  and  (-)  signs  correspond  to  the  left-  and  right-hand  sides  of  the 

branch  cut,  respectively.  Because  the  integrand  decays  as  —  on  the 

X. 

imaginary  axis,  it  is  not  possible  to  deform  the  contour.  Thus,  Eq.  (34) 
can  be  integrated  only  along  the  real  axis,  as  follows: 
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As  in  Eq.  (39),  the  first  integral  consists  of  both  real  and  imaginary  parts, 
and  is  integrated  within  a  finite  interval.  The  second  integral  has  only  a 
real  part  and  decays  at  least  as— as  X  -  co.  A  comparison  of  (43)  with 

A. 

Eq.  (39)  shows  that  the  second  integral  in  Eq.  (43)  converges  much  more 
slowly  than  that  in  Eq.  (39). 

/  j\ 

An  expression  for  f^;  can  also  be  obtained  (see  Appendix). 


III.  DIPOLE  ANTENNA  WITH  ONE -SIDE  EXCITATION 


As  shown  by  King  and  Wu  [8],  the  vector  potential  AZ(Z,R),  in  this 
case,  can  be  written  as  follows: 

1  +H  i 

AZ(Z,  R)  =  4^-  /  It(Z-)  K(Z-Z«;  R)  dZ<  -  -  V  Dq(Z;  R),  (44) 

o  -H  o-o 

where 


.  A  +7T 

D  (Z;  R)  =  /  R'  dR'  /  d0  - 


-i  V  Z2+(R2+R,2-2RR'  cos  0) 


0 


-7f 


a/ 


(45) 


Z2+(R,2+R2-2RR'co3  0) 


The  corresponding  integral  equation  for  the  total  current  distribution  is 


+H 

/  It(Z')  K(Z-Z';A)  dZ'  =—  (ccosZH-  |  sin  |z|)  +— VDq(Z;A); 
"H 


Z  ^H, 


(46) 


which  is  equivalent  to  Eq.  (la).  Since  now  the  5 -function  generator 
exists  only  on  the  outside  surface  of  the  antenna,  the  singular  current 
(which  is  later  subtracted  out  in  the  numerical  solution)  should  be  differ¬ 
ent  from  that  in  the  previous  case.  Aside  from  this,  Eq.  (44)  is  essen¬ 
tially  the  same  as  Eq.  (4)  except  that  it  is  necessary  to  evaluate  an 
extra  term  which  appears  in  the  column  matrix  ^  in  Eq.  (24a). 


A.  CURRENT  BEHAVIOR  NEAR  THE  DRIVEN  POINT 
For  an  infinitely  long  antenna,  i.  e.  ,  H  —  oo,  Eq.  (46)  degenerates  into 

f  I,(Z')K(Z-Z';A)dZ'  =  —  Ve+1lZl  +  — VD  (Z;A);  -oo<Z<oo.  (47) 

J  X  Y  r  O 

-CO 


3-1 


3-2 


Now,  define  the  Fourier  transform  of  Dq(Z;  R)  as 

D  (Z;  R)  =  f  D  (\;R)  e~lXZ  d\  (48a) 

o  J  o 

C 

1  /<+00  4-'\  7 

Do(X;  R>  =  j -  /  Do(Z;  R)  elX^dZ 

“CO 

=  fQA  R'  dR'  (t)  ^(k^2)  J0(R.V^) 

= - p -•  J^A's/  l-^2)  H^f  W  l-X2)  (48b) 

2  V  l-\2 

In  the  last  step  in  the  derivation  of  Eq.  (48b),  use  has  been  made  of  the 
relationship, 

Z 

/  tJ  (t)  dt  =  ZJJZ). 

Jo  0  1 

The  substitution  of  Eqs.  (5a-b),  (6a-b),  (7a-c)  and  (48a-b)  in  (47)  gives 

2*Mt)  h<i>(aV7^)  jo(aV7?) 


(49) 
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Now  apply  the  Wronskian, 

H^Z)  J Q(Z)  -  H ^(Z)  JX(Z)  =  ^ 
or 

i  w ,  M1}<z>  .iM 

H^fZ)  Jo(Z)  \H^'(Z)  Jo(Zy 


and  the  inverse  Fourier  transform  to  Eq.  (49), 

1  h(11)(a\/iT2) 


UZ)  =i^v  / 


O  C  \/ 1  -X2  h'^aVTT2) 


-i\Z 
e  dx 


(50) 


For  Z  >  0,  and  as  X  —  co  in  the  lower  half  plane. 


i  hJ^aVh  )  A 
A J 1  -x2  H^(aV  l-x2) 


1  -  I  X I Z 

e  1  ~  const  —  •  e  m 


Therefore,  the  contour  can  be  deformed  onto  both  sides  of  the  branch  cut, 


dX 


(51) 
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Now  apply  the  Wronskian  K^(Z)  H(2)(Z)  -  H(1)(Z)  H(2)(Z)  =-~to  Eq.  (51). 

10  0  1  7TZ/ 

The  result  is 


■i4 


iA  f  r°  r~ico\  e"iXZ 

I,  (Z)  =  —  V  (  J  +  /  '  ' 

1  V-1  J0 


'  A|V 


X2-l 


I  vs 


J*(AlV  ^|)  +  y2(a|V^7|) 


4V 


/ 


+iXZ 


dX 


'^0  l  0  t2 


J^(aVi-X2)  +  Y2(aV  1-X2) 


00 

-if 


-xz 


dX 


0  J2 
o 


(a\/i+x2)  +  Y2(aV  1+X2) 


(52) 


The  behavior  of  these  two  integrals  near  Z  =  0  has  been  discussed  by 
Duncan  [ll].  The  leading  term,  however,  can  easily  be  obtained  by  com¬ 
paring  Eq.  (50)  with  Eq.  (10).  Since  the  asymptotic  expressions  of  the 
two  integrands  in  Eq.  (50)  and  Eq.  (10)  behave  exactly  in  the  same  man¬ 
ner  as  X  —  co  on  the  real  axis,  the  leading  term  of  I^(Z)  in  Eq.  (50)  will 
then  be 


I  (Z)  ~  In  X  Z 

l  r  o 


as  Z  —  0. 


(53) 


Furthermore,  we  can  actually  prove  that  this  logarithmic  current  exists 
only  on  the  outer  surface.  That  is,  if  we  substitute  Eqs.  (6a,  b),  (48a,  b) 
into  Eq.  (44),  as  H  —  oo, 

Az(z)  =  if-  /  \M  h11)(r>aA:^)  J><V^K*Z 


+ 


A 


4v£, 


/  JjIaVh5)  h^ra/TY2) 


e"lXZdX  (54) 


o  C  Vl-X2 
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where  R>  =  {r,  A}.  A  comparison  of  (54)  with  Eq.  (25),  after  the 

< 

introduction  of  Ln(\)  as  the  Fourier  transform  of  I.^(Z),  leads  to 


I.  (X)  =  zillA  j 
in'  2 


in 


lt(\)  Vl  -X2  H^aVi?)  jo(aa/  i-x2) 


7TA' 

2^ 


VH^A  V l-X2)  Jx(aV  l-X2). 


(55) 


The  substitution  of  Eq.  (50)  into  Eq.  (55)  yields 


I.  (X)  s  0. 
in 


Therefore,  I-n(Z)  vanishes  everywhere  in  the  case  of  the  infinitely  long 
antenna,  and  Iq^(Z)  increases  logarithmically  near  the  origin. 

B.  NUMERICAL  SOLUTION  FOR  THE  CURRENT  DISTRIBUTION 


Let  D^.  =  D^(Z^.;  A).  With  Eq.  (43),  Eq.  (46)  leads  to 


o  K 


2N+1 


VUiaZ K +DK>  '  1  fK  ’  K=  1,  2, . . .  2N+1, 

m=l  o 


(56) 


where  and  fj^  are  defined  in  Eqs.  (17)  and  (23),  respectively. 

In  order  to  evaluate  that  part  of  the  current  which  is  distributed  on 
the  outside  (inside)  surface,  an  expression  for  the  difference  current  can 
also  be  obtained  by  the  previous  method.  However,  in  this  case,  the 
difference  current  also  has  a  logarithmic  singularity  at  the  driving  point. 

By  comparison  with  Eq.  (25),  the  magnetic  field  now  becomes 

-k  a  +H  ik  v  a 

B0(Z;  R>  *  8R  /  It(Z'>  K(Z-Z';R)  dZ'  +  gg  D0(Z;R)  (57) 
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Therefore, 

Id(Z)  =  27Tavo[B9(Z;  R=A+e)  +  BQ(Z;  R=A-e)] 
+H 


^  /  It(Z') 

2  -H  1 


K(Z-Z';R) 


R=A+e 


+  al  K(Z-Z';R) 


R=A-e 


dZ* 


i7TAV 


'o  i- 


9R 


R=A+e 


+  9R  Do(z;R) 


R=A-e 


as  e  -  0  and  Z  ^  H 


(58) 


With  the  definition. 


g(Z)  = 


i7TAV 


8R  Do(Z;R> 


R=A+e 


+  SR  Do(Z;E) 


-  —In 

|z|] 

R=A-e_ 

7T 

H  J 

as  e  -  0  and  Z  <  H 


(59) 


and  =  S(Zj^),  Eq.  (58)  can  be  transformed  into 
2N 

f,Kmm  2  AK  :  SK 


Ui.  1 

Tvw_  V  «(d)  T(t)xi  M)±  ~  ,  i2AV,  (K-l)A  T y  ,  „ 

It^  Lj  +  2  +  S^+  ..  In  — — ;  K  -  1 ,  2 , . . .  2N+ 1  (60) 


r(d) 

lK 


m=  1 


H 


where  pj^  and  f^P  are  defined  in  Eqs.  (28)  and  (29).  The  second  term 
in  Eq.  (59)  is  added  so  that  gK  is  finite  for  all  K. 


C.  EVALUATION  OF  DK  AND  gR 


From  Eqs.  (48a,  b), 
iA 


DK=if-  / - — J^aV l-\2)  H^aV l-\2)  e"lXKAd\ 

C  \/l-\2 


(61) 


3-7 


Aside  from  the  logarithmic  singularities  at  the  branch  points,  the  inte¬ 
grand  has  no  other  singularity.  Also,  as  |x|  —  oo  in  the  lower  half  plane. 


-iXKA 


Via2 


J^AV^1)  H^(aVI^)~  const -Le  ^ 


(62) 


so  that  the  contour  can  be  deformed  from  C  to  C1; 


dk  =  T  i(C+U 


0  /  e+i,r|Vx2-l| 


H^A  e+17r|V X2-l|)  dX 


0 


•ioc\ 


-iXKA 


+  I  /  +  /  ,  , _ 

1-1  0  ^  |V X2-l 


jx(a|a/  X2-l|) 


h^(a|a/ x2-i|)  dxl 


r  1  e+iKAX  /  j  o\  /  /  2\ 

iA  l  - _ r  J0(A/V  i-x2)  Jj(av  i-V) 


0  V^x2 


a  r~~—  J0(aV^)  j^an/I^5) 


e'KKAdX 


(63) 


0  v  l+x 


2  Z ! 

As  for  the  evaluation  of  gpr,  first,  define  L (Z)  =  ^ -*-n  ~ H-  ^or  1^1 
and  L(Z)  =  0  for  I  Z I  >  H.  Then  the  Fourier  transform  of  L(Z)  is 


+00 


LOO 


/  L(Z)  e+lXZ  dZ 


2ir 


■00 


H 


=  4  f  In  -§•  cos  XZ  dZ 

^  Jo  H 


S.(XH) 

7 rx 


(64a) 
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and 


HZ)  =  f  L(X)  e"lXZ  dX. 
C 


(64b) 


Here,  Si(XH)  is  the  sine  integral  whose  behavior  at  X.  —  oo  is  known.  It  is 


Si(XH)  -  -^cos  XH. 


(65) 


Now  substitute  Eqs.  (64a,  b),  (48a,  b)  into  Eq.  (59)  to  obtain: 


gK: 


iffAV 


f  HA) 


°  c 


k^an/Tv*)  +  -g-  ^ 


e-a(K-i)Adx 


(66) 


as  X  —  oo  on  the  real  axis. 


j^aV i-x2)  h^ax/  1-X2)  + 


i2 


Si(XH) 


7T2A  X 


•iXKA 


-iXKA 


77  AX 


+ 


77  A 


S-(i) *»(A) 


~  const 


-iX(K-l)A 


X 


(67) 


However,  since  Si(XH)  increases  exponentially  along  the  imaginary 
axis,  Eq.  (58)  is  not  true  on  both  sides  of  the  branch  cut,  and 
the  integration  has  to  be  performed  on  the  real  axis. 
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gK 


277 A2 vj  r1 

’  I  0 


J^aV l-\2)  Hf)(AV^)  + 


i2 


Si(XH) 


tt2A  X 


cosX(K-l)A  dX. 


—  f° 

77  J, 


Si(XH) 


77A  X. 


cos  X(K-1)A  dX. 


J 


(68) 


2 

Here,  g^  is  indeed  of  the  order  of  A  .  Therefore,  aside  from  the  driving 
point,  in  Eq.  (60)  is  equivalent  to  in  Eq.  (30),  in  the  case  of  thin 
antennas . 


IV.  NUMERICAL  RESULTS 


Based  upon  Eqs.  (24),  (27),  (56),  and  (58),  a  numerical  solution 
for  the  current  distributions  along  antennas  with  various  radii  and  heights 
can  be  obtained  for  both  1 -sided  and  2 -sided  excitation.  Also,  the  input 
conductance  of  the  antenna  is  readily  known  as  the  real  part  of  the  total 
current  at  the  driving  point,  in  the  case  of  2-sided  excitation  and  that 
of  the  outside  current  in  the  case  of  1 -sided  excitation.  The  input 
susceptance  of  the  antenna,  however,  is  always  infinite  in  either  case, 
due  to  the  fact  that  the  singular  current  near  the  driving  point  is  purely 
imaginary. 

Also,  in  the  process  of  numerical  computation,  the  choice  of  the 
number  of  sample  points  depends  upon  the  accuracy  required  for  the  solu¬ 
tion  on  the  one  hand,  and  the  computation  lime  on  the  other.  As  a  com¬ 
promise,  a  less  elegant  method  is  adopted.  That  is,  the  antenna  height 
is  first  divided  into  N  segments  (thus,  2N+1  sample  points),  and  a  set  of 
solutions  is  obtained  at  these  sample  points.  Then,  the  number  of  seg¬ 
ments  is  doubled  and  another  sec  is  obtained.  N  is  required  to  be  large 
enough  so  that  the  difference  between  these  two  sets  is  less  than  1-2%. 

In  this  mannc  although  it  is  not  known  exactly  how  this  approximate 
solution  converges  to  the  true  solution,  at  least  a  self-consistent  one  is 
obtained. 

Accordingly,  numerical  data  for  antennas  with  radii  ranging  from 

cl  h 

=  0.03907  to  0.  5,  and  antenna  heights  ranging  from-^-=  0.  1  to  0.  5  have 

cl 

been  obtained.  For  antennas  with  radii  ^-  <  0.  16,  the  input  conductances 
vs.  antenna  heights  are  compared  with  the  experimental  data  obtained 
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by  Holly  [12].  For  antennas  with  larger  radii,  no  experimental  data  are 
yet  available. 


A.  CURRENT  DISTRIBUTION 
1.  Half- Wavelength  Antenna  with  2 -Sided  Excitation 

Current  distributions  for  X- antennas  with  2 -sided  excitation  are 

shown  in  Figs.  (4a)-(4i).  The  real  part  of  the  total  current  I-p(Z)  be¬ 
haves  very  much  like  a  combination  of  terms  made  up  of  cosines  and 
shifted  cosines  with  half-angle  arguments,  except  near  the  ends  of  the 
antenna  where  the  slope  is  larger  than  given  by  these  two  terms.  This 
sharp  increase  in  slope  can  be  interpreted  as  a  consequence  of  an  ex¬ 
tremely  large  accumulation  of  charge  at  the  ends  of  the  antenna.  Further¬ 
more,  the  accumulated  charges  are  distributed  on  both  inside  and  outside 
surfaces  of  the  tubular  antenna;  they  are  proportional  to  the  slopes  of 
the  real  parts  of  the  outside  and  inside  currents,  respectively.  The  imag¬ 
inary  part  of  the  total  current,  on  the  other  hand,  can  be  approximated 
by  trigonometric  functions  only  after  the  subtraction  of  a  logarithmic 
term.  Also,  the  slope  of  this  part  at  the  ends  is  less  than  that  of  the 
real  part  in  most  cases.  It  follows  that  the  predominant  part  of  the 
charge  accumulation  is  proportional  to  the  time  rate  of  change  of  the  real 
part  of  the  current,  with  which  it  is  in  phase  quadrature. 

For  an  antenna  with  a  radius  less  than  the  cutoff  wavelength  of  the 

lowest  waveguide  mode  (i.  e. ,  X™.  ),  the  real  part  of  the  outside  cur- 

V01 

rent  can  be  approximated  by  a  combination  of  cosine  and  shifted  half¬ 
cosine  terms  with  an  argument  of  k^z+e),  where  e  is  a  certain  number 
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which  has  the  physical  interpretation  of  an  "end  correction."  At  the  end 
of  the  antenna,  this  real  part  of  IQut(Z)  remains  finite  and  reverses  its 
direction  as  the  current  passes  around  the  edge  of  the  cylinder  from  the 
outside  surface  to  the  inside  surface  of  the  antenna.  Once  inside,  it 
decreases  toward  the  driven  point.  The  rate  of  this  decrease  is,  in  gen¬ 
eral,  comparable  to,  but  certainly  not  the  same  as  that  of  the  TMq^ 
mode,  since  other  higher  modes  also  exist  near  the  end.  As  indicated 

in  the  figures,  a  is  the  "skin-depth"  of  the  TMq^  mode  where  the  cor- 

1  £1 

responding  current  decays  —  of  its  value  at  the  end.  As  becomes 
larger,  the  attenuation  decreases.  Also,  from  the  viewpoint  of  the  in¬ 
side  generator,  as  the  radius  of  the  "waveguide"  (i.  e.  ,  the  inside  of  the 
antenna)  becomes  larger,  radiation  due  to  the  open  ends  increases.  Thus 

3 

at  y-  =  0.  38  which  is  barely  under  the  cutoff  wavelength,  the  inside 
generator  actually  contributes  to  the  radiated  power,  as  the  real  part  of 
the  inside  current  becomes  positive  at  the  driving  points. 

On  the  other  hand,  for  half-wave  antennas  with  smaller  radius 
~  0.  04^  the  imaginary  part  of  the  outside  current  is  inductive  along 
most  of  the  antenna,  except  in  the  region  near  the  driving  point  where 
it  is  capacitive  with  a  logarithmic  singularity.  At  the  end  of  the  antenna, 
it  has  a  finite  inductive  value  on  the  outside  surface.  As  it  changes  direc¬ 
tion  it  becomes  a  capacitive  current  on  the  inside  surface,  where  it 
decreases  very  rapidly  toward  the  driving  point,  where  the  capacitive 
logarithmic  current  becomes  dominant  again.  Thus,  the  total  current 
distribution  along  half-wave  dipoles  that  are  not  too  thick  is  mostly 
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inductive.  As  becomes  larger,  the  range  over  which  the  imaginary 

part  of  the  total  current  is  inductive  is  reduced;  at  -^  =  0.  38,  it  becomes 

totally  capacitive.  At  the  same  time,  the  region  on  the  inner  surface  of 

the  tube  in  which  the  imaginary  part  of  the  current  decreases  toward  the 

driving  point  disappears  completely. 

For  an  antenna  with  a  radius  larger  than  \TM  ,  the  inside  gener¬ 
al 

ator  begins  to  play  an  active  role.  The  real  part  of  the  current  on  the 
inside  surface,  instead  of  having  a  maximum  at  the  end  of  the  antenna 
and  decreasing  toward  the  driving  point,  now  increases  toward  the  gener¬ 
ator  and  reaches  the  maximum  at  the  driving  point.  The  imaginary  part 
of  the  current  on  the  inner  surface  of  the  antenna,  on  the  other  hand,  is 
capacitive  when  ^-  =  0.45  and  inductive  when  —  =  0.  50.  As  will  be  shown 
later  in  section  C,  the  half- wavelength  antenna  is  not  yet  resonant  v/hen 
=  0.  45,  but  it  has  passed  through  resonance  when  =  0.  50.  In  either 
case,  the  magnitudes  of  both  the  real  and  imaginary  parts  of  the  current 
on  the  inner  surface  of  the  antenna  are  larger  than  their  counterparts  on 
the  outside. 

2.  Half-Wavelength  Antenna  with  One-sided  Excitation 

Current  distributions  for  X -antennas  with  1 -sided  excitation  are 

shown  in  Figs.  (5a)-(5i).  For  antennas  with  radii  less  than  XTM  ,  the 

01 

real  part  of  both  inside  and  outside  currents  behave  very  much  the  same 
as  with  2-sided  excitation.  At  —  =  0.  38,  the  real  part  of  the  current 

❖ 

Inductive  in  the  sense  that  the  region  where  the  imaginary  current  is 
inductive  is  much  larger  than  the  region  which  is  capacitive.  Of  course, 
in  any  case,  the  input  susceptance  is  always  capacitive. 
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remains  negative  everywhere  on  the  inside  surface  of  the  antenna,  instead 
of  becoming  positive  near  the  driving  point  as  it  does  with  2-sided  excita¬ 
tion.  This  reveals  that  the  inside  generator  actively  delivers  power  even 
though  the  interior  of  the  tube  is  still  well  under  the  cutoff  of  the  lowest 
waveguide  mode. 

The  imaginary  part  of  the  current  with  1 -sided  excitation  is  notice¬ 
ably  different  from  that  with  a  2 -sided  generator.  Specifically,  on  the 
inner  surface,  this  current  is  finite  everywhere.  There  is,  of  course, 
no  discontinuity  due  to  a  generator  ai  Z  =  0  on  the  inside  surface. 

For  an  antenna  with  a  radius  larger  than  ,  the  magnitudes 

M01 

of  both  the  real  and  imaginary  parts  of  the  current  on  the  inside  of  the 
antenna  increase  toward  the  generator  and  reach  a  maximum  at  the  driving 
point.  Also,  the  phase  relationship  between  the  real  and  imaginary  parts 
of  the  current  remains  constant  along  the  antenna  (except  in  the  region 
near  the  end),  since  they  both  form  the  standing-wave  pattern  character¬ 
istic  of  a  short-circuited  waveguide. 

3.  Full  Wavelength  Antenna  with  Two-sided  Excitation 

Current  distributions  for  the  full  wavelength-antenna  with  2-sided 
excitation  are  shown  in  Figs.  (6a)-(6g).  The  real  part  of  the  total  current, 
in  this  case,  can  still  be  approximated  by  a  shifted  cosine  and  a  shifted 
cosine  with  half-angle  argument,  while  the  imaginary  part  can  be  approx¬ 
imated  by  the  trigonometric  functions  only  after  the  subtraction  of  a  log¬ 
arithmic  term.  Near  the  end  of  the  antenna,  the  slope  of  the  real  part  of 
the  current  is  now  less  than  that  of  the  imaginary  part.  Therefore  the 
charge  accumulation  near  the  end  is  predominantly  real. 
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For  an  antenna  with  a  radius  less  than  XTM  ,  the  real  part  of  the 

iM01 

current  on  the  outer  surface  of  the  antenna  can  be  approximated  in  a 
fashion  similar  to  that  for  a  y  X- antenna.  However,  owing  to  the  length 
of  the  antenna,  this  real  part  of  the  current  now  is  negative  at  the  end. 

On  the  inner  surface,  the  real  part  of  the  current  decreases  rapidly  and 
there  exists  a  region  where  it  has  a  small  but  negative  value.  At  the 
driving  point,  however,  it  is  very  close  to  zero  and,  hence,  is  negligible 

Si 

compared  to  its  counterpart  on  the  outside  surface.  As  —  increases, 
this  region  with  a  small  negative  value  decreases  and  the  real  current  on 
the  inside  surface  becomes  more  and  more  positive.  At  —  =  0.  38,  it  has 
a  large  positive  value  everywhere  inside  the  antenna,  and  the  inside  gen¬ 
erator  again  contributes  significantly  to  the  radiated  power.  Also,  as 
y  increases,  the  imaginary  part  of  the  total  current  changes  from  mostly 
inductive  to  purely  capacitive,  while  the  current  on  the  outer  surface 
remains  inductive  and  is  not  sensitive  to  an  increase  in  y 

For  an  antenna  with  radius  larger  than  XTM  ,  the  inside  generator 

iiVi01 

again  plays  an  active  role.  The  real  part  of  the  current  on  the  inside 
surface  now  has  a  maximum  at  the  driving  point,  decreasing  toward  the 
end,  where  it  becomes  negative.  Also,  unlike  the  situation  in  the  y  X- 
antenna,  the  imaginary  part  of  the  current  on  the  inside  surface  is  mostly 
inductive  for  —  =  0.  45  and  0.  50.  This  is  a  consequence  of  the  fact  that 

K 

the  antenna  is  very  much  longer  than  resonant. 

4.  Full-Wavelength  Antenna  wim  One-sided  Excitation 

Current  distributions  for  the  X -antenna  with  1 -sided  excitation  are 
shown  in  Figs.  (7a)-(7g).  Similar  to  the  case  of  the  y  X -antenna,  the 


FIG  7  CURRENT  DISTRIBUTION  ON  FULL  -  WAVELENGTH  ANTENNA  WITH  l-SIDE  EXCITATION 
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distribution  of  the  real  part  of  the  current  behaves  very  much  like  that 


a  . 


with  the  1 -sided  excitation,  when  y  is  less  than  ,  and  the  imagin¬ 


ary  current  has  no  logarithmic  term  on  the  inside  surface.  For  y larger 

than  XTM  ,  away  from  the  end  of  the  antenna  both  real  and  imaginary 
iM01 

parts  of  the  current  on  the  inner  surface  follow  the  same  standing-wave 
pattern.  The  phase  relation  between  these  two  components  of  current 
also  remains  constant  along  the  antenna. 

B.  INPUT  CONDUCTANCE 


The  input  conductance  of  an  antenna  with  various  values  of  radius 
vs.  the  antenna  height  are  shown  in  Figs.  (8a) -(8c).  For  convenience 
in  the  description,  let  the  whole  antenna  be  treated  as  consisting  of  two 
coupled  radiating  systems:  system  (T)  is  driven  by  the  generator  on  the 
outside  surface;  system  (2)  by  the  generator  on  the  inside.  The  two  sys¬ 
tems  are  then  coupled  by  way  of  the  open  end  of  the  antenna.  Thus,  cur¬ 
rents  and  voltages  at  the  two  driving  points  can  be  related  by  the  equation 


h-YnYi  +  Yi2Y2 


In  =  Yr,,  V.  +  YnoV 


(69) 


21  v  1 


22  2 


where  and  Y are  the  self-admittances  of  systems  (T)  and  (2) ,  re¬ 
spectively,  and  Y-^2  and  Yg^  are  the  mutual  admittances.  The  input  con¬ 
ductances  G1  and  G2  of  an  antenna  with  1 -sided  and  2-sided  excitation 
are  given,  respectively,  by 


G(1)  =  Real  part  of  Yn  =  Gn 

i(2) 


LG  =  Real  Part  0f  (YirY22+Y12+Y21)  =  Gll  +  G22  +  G12  +  G 


(70) 

22 
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(1)  Antenna  with  moderately  large  radius  <  0.  16  or  277a  <  X, 
i.  e. ,  the  circumference  of  the  antenna  is  less  than  a 

For  antennas  in  this  range,  the  radius  is  much  smaller  than  the 
cutoff  wavelength  of  the  TMq^  waveguide  mode,  even  though  it  may  be 
very  large  compared  to  the  conventional  dipole  antenna.  Therefore,  the 
input  conductance  is  not  sensitive  to  the  method  of  excitation  (i.e.,  » 

G12,G22^*  The  general  behavior  of  the  input  conductance  vs.  antenna 

cl 

height  is  also  similar  to  that  of  the  conventional  thin  antenna.  As  —  is 

increased,  both  values  of  the  input  conductance  at  resonance  (G_^  ) 

r  max 

and  antiresonance  (G  .  )  increase,  but  the  ratio  of  the  two  decreases. 

mm 

Also,  the  resonant  and  the  antiresonant  lengths  decrease.  Thus,  the 
corresponding  curve  in  Fig.  (8a)  will  move  up,  shift  toward  the  left  and 
become  flatter. 

Data  obtained  from  the  numerical  solution  are  in  good  agreement 
with  those  obtained  experimentally  by  Mack  [  1 3] ,  in  the  case  of  conven¬ 
tional  thin  antennas.  For  antennas  with  moderate  large  radii  0.03907^ 
the  agreement  with  data  obtained  by  Holly  [  1 2]  is  also  reasonably  good. 

The  two  curves  seem  to  differ  by  less  than  0.  0075  X,  and  the  experimental 
value  at  resonance  is  consistently  less  than  the  numerical  solution  by 
only  3-4%.  (For =  0.01074,  the  comparison  is  actually  based  on  an 

3.  cl 

average  experimental  value  of  —  =  0.  10157  and  — =  0.  11287,  therefore 
is  less  accurate.) 

(2)  Antenna  with  radius  less  than  X-p^.  (0.  16  <  -—  <  0.  383): 

As  the  radius  of  the  antenna  becomes  larger,  the  difference  between 


'  Resonance  (or  antiresonance)  is  defined  as  the  position  where  the  input 
conductance  reaches  maximum  (or  minimum). 
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the  input  conductance  with  1- sided  and  2 -sided  excitations  also  becomes 

more  significant.  For  a  fixed  radius,  the  corresponding  curve  for 

(2) 

becomes  sharper  than  that  of  G  .  Also,  it  has  a  larger  resonant  con¬ 
ductance,  but  smaller  antiresonant  conductance.  The  resonant  lengths 
for  the  two  different  excitations  are  different,  too,  since  the  behavior  of 
G^  is  no  longer  approximately  the  same  as  (G^+G^'^-l+C^^4  *Purt*ier" 
more,  when  the  antenna  is  short,  the  coupling  between  the  inside  and  the 
outside  surfaces  is  comparatively  stronger,  and  the  real  part  of  the  cur¬ 
rent  on  the  inner  surface  is  negative  at  the  driven  point.  When  the  an¬ 
tenna  is  long,  the  coupling  is  comparatively  weaker,  ana  the  real  part  of 
the  current  on  the  inner  surface  is  positive  at  the  driven  point.  There¬ 
fore,  G^  is  quite  a  bit  larger  than  G^  when  the  antenna  is  short,  and 

(2) 

is  a  little  smaller  than  G  when  the  antenna  is  long. 

As  the  radius  of  the  antenna  is  increased,  both  the  input  conduct¬ 
ances  at  resonance  G^  (or  G^  )  and  at  antiresonance  „  (or  G^|  ) 

max  max  mm  mm 

become  larger.  However,  the  resonant  lengths  now  increase  and  the 


antiresonant  lengths  remain  relatively  unchanged. 

The  behavior  of  G^  and  G^  near  resonance  and  antiresonance, 


and  also  the  location  of  resonance  as  a  function  of  y  are  shown  in  Fig.  9. 

(3)  Antenna  with  radius  larger  than  >  0*  383^: 

Since  the  TMq^  mode  can  now  propagate  in  the  inside  of  the  antenna, 

both  G^  and  G^  reach  a  sharp  maximum  when  this  waveguide  mode  is 

(2) 

at  resonance.  However,  for  a  fixed  radius,  the  value  of  G'  is  much 

llla-X 

larger  than  that  of  G^  .  This  is  due  to  the  fact  that  the  TMni  mode  is 

max  0 1 

driven  at  the  location  of  the  maximum  current  with  2 -sided  excitation. 


but  excited  only  by  the  coupling  from  the  outside  at  the  end  of  the  antenna 
with  1 -sided  excitation. 
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FIG.  9.  INPUT  CONDUCTANCE  AND  ANTENNA  HEIGHT  AT  RESONANCE 


V.  CONCLUSION 


The  problem  of  the  thick  dipole  antenna  has  been  investigated 
theoretically  for  both  1-sided  and  2-sided  excitations.  When  the  circum¬ 
ference  of  the  antenna  is  less  than  one  wavelength,  the  solution  obtained 
from  the  physically  more  realizable  model  with  1- sided  excitation  does 
not  differ  significantly  from  the  conventional  model  with  2-sided  exci¬ 
tation,  However,  for  antennas  with  larger  radius,  it  is  only  correct  to 
use  the  model  with  1 -sided  excitation  in  order  to  obtain  a  solution  that 
is  comparable  with  measured  results. 
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APPENDIX 


Define 


+2A  .1  , 

P(m,X)  =  X  /  In  \~\  e"1lmA-Z'|\  dZ, 

-2A  ^a 

2A 

P(0A)  =  2\  /  ln^e“iZ,XdZ' 

=  ^^(e-12^-!)  ^  -  /2A  dZ,' 

L  o  o  ^ 

=  -2i[Ci(2XA)  -  v  -  In  2XA+  iSi(2XA)] 

P(1A)  =  X  I  ln  H&I  e'l(A"Z')X  dZ'  +  X  f 2 ^ in ~~  e*i(A-Z')X 


(A-l) 


=  2X  cos  XA  f  ln^  e~lZ  X  dX  +  i2X  /  In  ^  sin  (Z'-A)X  dZ' 


=  [Ci(2XA)  -  v  -  In  2XA  -iSi(2XA)] 


+  2iK  -  J In  ~  dtcosKZ1  -  1)  -  imAX  f  la& 


d  sin  Z'X 


j.  ln  2  ,  A 

+  — r—  cos  (Z'-A)X 

0 


2i  {cos  XA[  -Ci(2XA)  +  Ci(XA)  +  iSi(2XA)]  +  sinXASi(AX)  +  In  2}  (A- 2) 
2a 

P(m,X>  =  \  /  ln-g  [e-i<mA-Z')X  +  e-i(mA+Z')\]  dz, 


-imAX 


r  Z1 

J  In  ^  d  sin  Z'X 


_2e-im^X  Si(2AX). 


m  >  2, 


(A-3) 


A-l 


A-2 


where  Ci(Z),  Si(Z)  are  the  cosine  and  sine  integrals,  respectively.  Sub¬ 
stitute  Eqs,  (A-l,  -2,  -3)  into  Eqs.  (17),  (29)  and  compare  with  Eqs.  (39) 
and  (43): 


v  in  (A-l)  is  the  Euler's  constant. 
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